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We investigate the critical properties of the 5 = 1/2 bond-alternating spin chain with a next- 
nearest-neighbor interaction in magnetic fields. By the numerical calculation and the exact solution 
based on the effective Hamiltonian, we show that there is a parameter region where the longitudinal 
incommensurate spin correlation becomes dominant around the half-magnetization of the saturation. 
Possible interpretations of our results are discussed. We next investigate the effects of the interchain 
interaction (J'). The staggered susceptibility and the uniform magnetization are calculated by 
combining the density-matrix renormalization group method with the interchain mean-field theory. 
For the parameters where the dominant longitudinal incommensurate spin correlation appears in 
the case J' = 0, the staggered long-range order does not emerge up to a certain critical value of J' 
around the half-magnetization of the saturation. We calculate the static structure factor in such a 
parameter region. The size dependence of the static structure factor at = 2kp implies that the 
system has a tendency to form an incommensurate long-range order around the half-magnetization of 
the saturation. We discuss the recent experimental results for the NMR relaxation rate in magnetic 
fields performed for pentafluorophenyl nitronyl nitroxide. 

PACS numbers: 75.10.Jm; 75.40.Cx; 71.10.Pm; 76.60.-k 



I. INTRODUCTION 

Onc-dimcnsional (ID) spin-gapped systems have at- 
tracted a great amount of attention both theoretically 
and experimentally. It was shown for the ID S = 1/2 
spin-gapped systems that a noticeable feature of each 
system appears in critical properties, when the energy 
gap is collapsed by external magnetic fields 0, In 
ffci < H < , where and are the lower criti- 
cal field and the saturation field respectively, the field de- 
pendence of the critical exponent of the spin correlation 
function was calculated numerically for the Haldane-gap 
system 0, Q , the S = 1/2 bond-alternating chain , the 
alternating-spin chain [g , and the S = 1/2 two-leg ladder 
0, 13 ■ It was shown that in these systems the transverse 
staggered spin correlation makes a leading contribution, 
and that the critical exponent exhibits characteristic be- 
havior as function of magnetic field or the magnetiza- 
tion in each system. Such features can be observed by 
the field dependence of the divergence exponent of the 
NMR relaxation rate l/Ti with decreasing temperature 
Experimentally, the divergence exponent of 1/Ti 
was measured in the 5 = 1/2 two-leg ladder and the 
Haldane-gap system The results were discussed in 
connection with the theoretical results Q, 0, S S Q • 

When temperature is further decreased and the in- 
terchain or interladder interaction becomes relevant, a 
quantum phase transition towards a three-dimensional 
(3D) ordered state takes place. Such 3D ordered states 
in magnetic fields were observed experimentally in sev- 
eral quasi-lD spin-gapped systems such as Haldane-gap 
materials NDMAZ O and NDMAP [pL a 5 = 1/2 
two-leg spin-ladder material CuHpCl [I^, and S = 
1/2 bond-alternating chains Cu(N03)2 ■ 2.5H2O ^ and 
(CH3)2CHNH3CuCl3 [13. 



The S = 1/2 bond-alternating spin chain with a 
next-nearest-neighbor (NNN) interaction is also a typi- 
cal ID spin-gapped system. Fascinating phenomena have 
been investigated intensively using this model. It was 
shown numerically that there appears a plateau region on 
the magnetization curve at half of the saturation value 
|l8l Il9j . The phase diagram was determined precisely 
by the level spectroscopy analysis Using bosoniza- 
tion technique and a perturbation calculation, a simple 
picture of the half-magnetization-plateau state is pre- 
sented as the twofold degenerate state with the singlet 
and triplet pairs occupying the strong bonds alternately 

M- 

This model can be regarded as a minimal model for 
the spin-Peierls material CuGeOs |2lLl23. Note that the 
parameter sets proposed for CuGeO'^ |2ll |2^ lie within 
the half-magnetization-plateau region of the phase dia- 
gram 19]. It was shown for these parameters that in 

< H the static spin susceptibility parallel to mag- 
netic fields takes the maximum at the incommensurate 
(IC) wave vector Using the adequate parameters 
for CuGe03 [22, the critical exponents of the spin cor- 
relation functions were further investigated numerically 
in iJci < H < 24]. It was shown that the leading 
contribution of the spin correlation function changes de- 
pending on magnetic fields: In the middle range between 

< H < the IC spin correlation parallel to mag- 
netic fields becomes dominant, while around H = Hc^ 
and H = Hc2 the staggered spin correlation perpendic- 
ular to magnetic fields becomes dominant. Such critical 
properties may have a relation to the IC phase of CuGe03 
in magnetic fields Furthermore, the results suggest 
that the IC long-range order may be stabilized, when the 
interchain interaction is relevant. 

In this paper, we investigate critical properties of a 
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S = 1/2 bond- alternating spin chain with a NNN in- 
teraction in magnetic fields. We turn our attention 
to the dominant longitudinal IC spin correlation. We 
then investigate the characteristics of the ordered state, 
when the interchain interaction is taken into account. 
In Sec. II, we calculate critical exponents of the spin 
correlation functions systematically in various parame- 
ter sets, combining a numerical diagonalization method 
with finite-size-scaling analysis based on conformal field 
theory |^. On the basis of the results, we determine 
the phase diagram for the dominant longitudinal IC spin 
correlation and discuss the origin of such properties. In 
Sec. Ill, we next investigate the ordered state in mag- 
netic fields, by taking account of the interchain inter- 
action. Combining the density-matrix renormalization 
group (DM RG) method with the interchain mean-field 

theory [H m El m US m 113 , we calculate the stag- 
gered susceptibility and the uniform magnetization. We 
discuss whether the staggered long-range order can be 
stabilized in < H < Hc^. To investigate characteris- 
tics of the ordered state, we calculate the static structure 
factor. In Sec. IV, we discuss the recent experimental re- 
sults for 1/Ti in magnetic fields performed for pentaflu- 
orophenyl nitronyl nitroxide (F5PNN) Sec. V is 

devoted to the summary. 

II. CRITICAL PROPERTIES 

A. Model and method 

Let us first consider the ID S = 1/2 bond-alternating 
spin system with a NNN interaction in magnetic fields 
described by the following Hamiltonian, 

n^Ho + Hz, (1) 



Hamiltonian into the subspace according to the wave vec- 
tor k and the magnetization M — J2i=i + ^i,r^ ■ 
The distance between the neighboring two unit cells is 
set to unity. Thus, the wave vector takes the discrete 
value k = {2tt/L) x integer for finite L. Using Lanczos 
algorithm, the lowest energy in each subspace is calcu- 
lated numerically. For the 2L-spin system, we define the 
lowest energy of in the magnetization M and the wave 
vector k as {L, M). In given L and M, Ek {L, M) takes 
the minimum at A; = feg. We simply describe Ekg{L,M) 
as E{L,M). Following the method developed in Ref. 3, 
we investigate critical properties of the system in mag- 
netic fields. 



B. Central charge 

We first investigate the central charge c by use of 

li?(i,M)^.(m)-^-l, (4) 

where e(m) is the ground state energy per a unit cell in 
the thermodynamic limit with m — M/ {2L) (0 < m < 
1/2). The velocity is estimated as 

Vs = ^[EkAL,M)-E{L,M)], (5) 

where ki is the wave vector closest to feg . From the L de- 
pendence of E{L, M)/ L, we derive (7r/6)CTs. Combining 
(tt/6)cVs thus obtained and numerically calculated Vs, we 
can evaluate the central charge c numerically. Typical re- 
sults are shown in Table I. From the results, we conclude 
that c = 1 for a = 0.10,0.15,0.18 and 0.20 in 6 = 0.15 
and the system in the gapless region can be described as 
the Tomonaga-Luttinger (TL) liquid. 



Ho = 2J^ Ul + (5) Su ■ S,,r + {l-S) S,,r ■ S.,+i,l 
i=l 

+ Oi[Si^l ■ Si-^.l^l + Si^r ■ Si+i^r] (2) 



L 

Hz = -9fiBHY,{Sli + Sl,), (3) 

2 = 1 

where L is the total number of unit cells, which consist of 
neighboring two spins, Si^K^r) is the 5* = 1/2 spin operator 
of the left- (right-) hand side in the ith unit cell, and 
H is the magnitude of the external magnetic field along 
the z axis. Here, S is the bond-alternation parameter 
with < S < 1 and 2aJ is the NNN antiferromagnetic 
coupling. We set J — gj-iB = 1- The periodic boundary 
condition is applied. 

Since the system has translational symmetry and ro- 
tational symmetry about the z axis, we can classify the 



TABLE I: The central charge for a = 0.10,0.15,0.18 and 

0.20 in (5 = 0.15^ 

a 0.10 0.15 0.18 0.20 

m = 1/8 1.10 1.10 1.10 1.12 

1/7 1.08 1.06 1.09 1.09 

1/6 1.07 1.06 1.08 1.08 

1/5 1.09 1.09 1.09 1.09 

3/14 1.06 1.06 1.06 1.06 

1/4 1.06 1.06 1.05 1.05 

2/7 1.05 1.06 1.05 1.05 

3/10 1.07 1.06 1.05 1.04 

1/3 1.06 1.06 1.05 1.05 

5/14 1.06 1.00 1.03 1.02 

3/8 1.07 1.06 1.05 1.03 

2/5 1.07 1.06 1.05 1.03 



It was shown that in the adequate parameters of the 
Hamiltonian Q the plateau region appears at half of the 
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saturation value (m = 1/4) in the magnetization curve 
[TsL l20| . We investigate the magnetization curve us- 
ing the asymptotic forms of the excitation energy 



E{L, M + 1) - E{L, M) - H+{m) - irv^ri' 



I 

T'' 



(6) 



E{L, M) - E{L, M-l)-H- (to) 



-7rw,77^-, (7) 



where rj^ is the critical exponent of the spin correlation 
function transverse to the magnetic field, H^{m) and 
H~{m) are the magnetic field in a given to in the ther- 
modynamic limit. If H~^{mo) — H~(mo), no plateau 
appears at to = toq in the magnetization curve. 

In the parameters used in Table I, H^{m) and H~{m) 
extrapolated into L — > cx) satisfy H^{m) = H^{m) 
within the accuracy of 0(10^"*). Therefore, we conclude 
that no plateau appears at to = 1/4 within our numerical 
accuracy. 



C. Critical exponents 

We next investigate the critical exponents of the spin 
correlation functions. The spin correlation functions of 
the long-distance behavior in the TL liquid take the forms 
as 



r^m^ + Cin-^ + C2n-'i' cospfcpn]. 



(8) 



(9) 



where cj)n ^'' is composed of the two spin operators at the 
nth unit cell and 2A:f — 7r(l — 2to). The critical exponents 
are given by 



L 



V 



[E2kAL,M)-E{L,M)]. 



(10) 



L 



[E{L, M + 1) + E{L, M - 1) - 2E{L, M)] . 

(11) 

Using the expressions (|10|l and (|ll|l . we calculate 77^ and 
T]^. Since the size dependence of 77^ and 77^ is found to 
be well fitted by 0{L~^), we extrapolate the results in a 
given m to L —> 00. 

In Fig. 1, we show the extrapolated 77^ and 77^ for 
a = 0.10,0.15,0.18 and 0.20 in 6 = 0.15 as function of 
TO. It is a difficult issue to obtain the precise results in 
TO < 0.1, because we cannot treat larger size systems. 
At the lower critical field and the saturation field, the 
system may be described by a boson with infinitely large 
repulsion. Therefore, the critical exponents tafte the val- 
ues 77^ = 2 and 77^ = 1/2 at to = and 1/2. On the basis 
of these results, we conclude that in a = 0.10, and 0.15, 
the relation 77^ < 1 < r;^ is satisfied in < to < 1/2. 



TABLE IL The value 7?^ • 77^ for a = 0.10, 0.15, 0.18 and 0.20 
in 5 = 0.15. 



a 0.10 0.15 0.18 0.20 

1/8 0.98 0.99 0.98 0.97 

1/7 0.99 0.99 0.99 0.99 

1/6 1.00 1.00 1.00 1.00 

1/5 1.00 1.00 0.99 0.99 

3/14 1.00 1.00 1.00 1.00 

1/4 0.99 0.96 0.93 0.91 

2/7 1.00 1.00 1.00 0.99 

3/10 1.00 0.99 0.99 0.99 

1/3 1.00 1.00 1.00 1.00 

5/14 0.99 1.00 1.00 0.99 

3/8 1.00 1.00 0.99 0.99 

2/5 1.00 1.00 1.00 0.99 
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FIG. 1: The extrapolated 77^ and 77^ for a = 0.10,0.15,0.18 
and 0.20 in S — 0.15 as function of m. 



The results indicate that the transverse staggered spin 
correlation is dominant in < 771 < 1/2. In a = 0.18 and 
0.20, on the other hand, there appears the region where 
the relation 77^ < 1 < 77^ is satisfied in 0.16 < to < 0.22 
and 0.13 < to < 0.28, respectively. The results indicate 
that the dominant spin correlation changes depending 
on the magnetization: For a = 0.18 the longitudinal IC 
spin correlation becomes dominant in 0.16 < to < 0.22, 
while the transverse staggered spin correlation becomes 
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dominant in < m < 0.16 and 0.23 < m < 1/2. For 
a = 0.20 the longitudinal IC spin correlation becomes 
dominant in 0.13 < m < 0.28, while the transverse stag- 
gered spin correlation becomes dominant in < < 0.13 
and 0.28 < m < 1/2. To check the numerical accuracy, 
we evaluate the value r]^ ■ rj^ . From the results shown in 
Table II, the universal relation r/^ ■ rj^ = 1 characteristic 
of the TL liquid 34, 35] is well satisfied in < m < 1/2 
except for m = 1/4 at a = 0.15, 0.18, and 0.20. 

The appearance of the half-magnetization plateau in 
this system is the Berezinskii-Kosterlitz-Thouless quan- 
tum phase transition |l9j . Therefore, in finite systems, rj^ 
and Ty^ at m = 1/4 are suffering from slowly-converging 
logarithmic size-corrections in the vicinity of the transi- 
tion point. As will be seen in Fig. 2, the parameters 
a = 0.15,0.18, and 0.20 in ^ = 0.15 lie close to the tran- 
sition point. It is thus difficult to obtain their highly 
accurate results for 77^ and 77^ at m = 1/4. 

We develop the calculation in other parameter sets of a 
and S, turning our attention to the behavior 77^ < 1 < 77^. 
The results are summarized in Fig. 2. In the darkly 
shaded area, the relation 77^ < 1 < 77^ is satisfied around 
the half-magnetization of the saturation. In the lightly 
shaded area, the relation 77^ < 1 < 7;^ is satisfied in 
< 771 < 1/2. The half-magnetization plateau appears 
in the 'plateau' area. The two broken lines are the bound- 
aries for the half-magnetization plateau obtained by the 
level spectroscopy analysis [l3|- Note that the method 
used in this paper is difficult to detect the small half- 
magnetization plateau close to the broken lines |36J. 

Such a dominant longitudinal IC spin correlation is 
closely related to the formation of the half-magnetization 
plateau, which can be regarded as the Akp CDW state 
|20|. It was shown that in the vicinity of the CDW 
transition point the dressed charge is suppressed and 
takes the value typical of the strongly correlated system 
mUllsiElill- In the ID Hubbard model, for exam- 
ple, the dressed charge of the charge excitation takes the 
value for the spinless fermion at half-filling irrespective 
of the strength of the Coulomb repulsion [s^, IH, Isof . 
Within the method used in this paper, it is difficult 
to investigate the critical properties close to the half- 
magnetization plateau. To see the critical properties in 
the other viewpoint, we calculate 7;^ in a ~ and (5 ~ 1 
on the basis of the effective model. 



D. Critical properties based on the effective model 

Under the condition a and 6^1, the Hamiltonian 
can be mapped onto the ID S —1/2 XX Z model in 
effective magnetic fields [H [l^ El : 

Woff' — Jcs l^'S'f S',f_^]^ + SfSf^i + AS^S^_^i 

i 

- iJeff$^5^ (12) 




FIG. 2; In the darkly shaded area, the relation 77^ < 1 < r;^ 
is satisfied around the half-magnetization of the saturation. 
In the lightly shaded area, the relation t;^ < 1 < 77^ is satis- 
fied in0<m<l/2. The half-magnetization plateau appears 
in the 'plateau' area. Two broken lines are the boundaries 
for the half-magnetization plateau obtained by the level spec- 
troscopy analysis The open circle corresponds to the 
parameters for F5PNN (a = 0.15 and 5 = 0.43). The four 
open triangles correspond to the parameters used in Fig. 1 
(a = 0.10, 0.15, 0.18 and 0.20 in 5 = 0.15). 

where Sf {a = x, y, z) is the pseudo spin operator made 
up of the two states of the two spins on the strong bond. 
The effective coupling constants and the effective mag- 
netic field are given by using the original parameters 
as Jeff = [2a - (1 - S)], AJ^s = [a + i(l - S)] and 
H,s = H-[a + ^{5 + 3S)]. 

When A > 1 with Joff > 0, the ground state has Neel 
order and the excitations are gapped. Since the mag- 
netization of the effective Hamiltonian 77icff satisfies the 
relation TT^eff -I- 1/2 = 2m, the ground state nics — cor- 
responds to the half-magnetization-plateau state of the 
Hamiltonian The critical exponents of the spin cor- 
relation functions in Tigw are the same as those in the 
Hamiltonian 0,E1|- Therefore, we calculate the 

critical exponents in TioS to see the critical properties 
around a ~ and 6 ^ 1 of the Hamiltonian 

The effective Hamiltonian enables us to calculate rj^ 
exactly in the gapless region using the Bethe ansatz so- 
lution E^ . The critical exponent nf is obtained from the 
dressed charge as 77^ = 2[2'(A)]^ [iy. The dressed charge 
Z{X) is obtained from the integral equation, 

1 

Z{\) = 1-— d\'a2{\-\')Z{\'), (13) 

where a„(A) = (/> sinh(7T,(/))/[cosh(r7,(/)) — cos((/)A)] with 
A = cosh (0 > 0) . The cutoff A is determined by the 
condition for the dressed energy £(±A) = 0, where e(A) 
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^eff 

FIG. 3: Heff dependence of rj^ and rries for several A. 



E. NMR relaxation rate 

We have shown that the longitudinal IC spin correla- 
tion becomes dominant around the half-magnetization- 
plateau region in the system described by the Hamilto- 
nian Q • Experimentally, such a feature can be observed 
by the NMR relaxation rate 1/Ti. When the NMR is 
done on the nuclei located at the different sites from the 
electronic spins, the relaxation occurs through a dipolar 
interaction between the nuclear and electronic spins. In 
this case, 1/Ti of the TL liquid is expressed as a sum 
of contributions from the longitudinal and transverse dy- 
namical spin susceptibilities. Paying attention to the di- 
vergence behavior of l/Ii, we obtain the expression 
l/Ti ^ C'2T''"-i + i)iT''"-\ where the first term orig- 
inates in the longitudinal IC spin susceptibility and the 
second one originates in the transverse staggered spin 
susceptibility. Note that C2 and Di depend on tempera- 
ture and magnetic field. However, their effects are prob- 
ably weaker than the divergence behavior. Therefore, we 
consider only the field dependence of the divergence ex- 
ponent of 1 /Ti . Depending on 77^ and rj^ , 1 /Ti takes the 
form as 



1 



(16) 



is obtained from the integral equation in given magnetic 
field as 

e(A) =eo(A) - — / dX'a2{X ~ X')e{X') (14) 

J -A 

with £o(A) = Heff — [27rJsinh(/)/0]ai(A). The magnetiza- 
tion is obtained by use of the dressed charge as 

1 1 

^^«=n-7r dXai{X)Z{X). (15) 

We calculate rj^ and nics as function of Hcs- The re- 
sults are shown in Fig. 3 for several anisotropic param- 
eters A. Around both ends of the half-magnetization 
plateau, the relation 77^ < 1 is satisfied and the lon- 
gitudinal IC spin correlation becomes dominant. It is 
calculated analytically that 77^ — 1/2 at the lower crit- 
ical fields i?eff,ci = ±[2 sinh{(/))/ (J)] K{m')m' and rj^ = 2 
at the saturation fields -ffeff,c2 — ±(1 + ^) irrespective of 
A(> 1), where K{m') is the complete elliptic integral and 
its modulus m! . The half-magnetization-plateau state in 
the Hamiltonian Q is Neel ordered in the language of the 
pseudo spin. As mentioned before, the dressed charge is 
suppressed in the vicinity of the CDW transition point. 
The results that 77^ takes the minimum at H^s — -f^cff,ci 
can be explained in this point of view. 

Judging from these findings, we conclude that such 
a dominant longitudinal IC spin correlation emerges 
around the half-magnetization-plateau region. 



We discuss the divergence property of 1/Ti using the 
results for 77^ and 77^ shown in Fig. ^ For a = 0.10 and 
0.15 in (5 = 0.15, the relation 77^ < 1 < 7;^ is satisfied in 
< 771 < 1/2, indicating that the divergence behavior is 
caused by the transverse staggered spin correlation. For 
a = 0.18 and 0.20 m 5 = 0.15, the relation 77^ < 1 < 77"^ 
is satisfied in 0.16 < m < 0.22 and 0.13 < m < 0.28, 
respectively. In the following, we express these regions as 
nip-^ < m < mp^. Thus, for a — 0.18 and 0.20 inS — 0.15, 
the relation 7;^ < 1 < 77^ is satisfied in < ttz < nip^ 
and uip^ < m < 1/2. Therefore, in these parameters the 
divergence behavior is caused by the longitudinal IC spin 
correlation in uip^ < m < nip^ , and by the transverse 
staggered spin correlation in < 771 < nip^ and 771^3 < 
m < 1/2. 

The m dependence of the divergent exponent 7 is sum- 
marized in Fig. 4. For a — 0.10 and 0.15 in 5 = 0.15, 7 
varies concavely as function of m with 7 = 1/2 at 771 = 
and 1/2. For a = 0.18 and 0.20 in S = 0.15, as 771 in- 
creases, 7 decreases in < 777 < tti^^ , 7 varies convexly in 
TTipj < 777 < TOj,2, and 7 increases in 77ip2 < 777 < 1/2. It is 
noteworthy that 7 = at 771 = 771^^ and TTip^, where 1/Ti 
shows no divergence and becomes almost independent 
of temperature. Such a feature can be observed experi- 
mentally. In fact, it was reported for F5PNN that 1/Ti 
exhibits power-law behavior in Hc^ < H < Hc^ , and at a 
certain magnetic field 1/Ti became almost independent 
of temperature 33|. The results may be an evidence for 
the change of the dominant spin correlation in magnetic 
fields. Details will be discussed in Sec. IV. 
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m 

FIG. 4: The divergence exponents 7 for a = 0.10,0.15,0.18 
and 0.20 in 5 = 0.15 as function of m, where 1/Ti ~ T~"' . 
The hnes are to guide the eyes. 



III. EFFECTS OF THE INTERCHAIN 
COUPLING 



A. Mean-fleld approximation 

We next consider the effects of an interchain interac- 
tion on the system described by the Hamihonian 
The Hamiltonian may be written down as 

^ = { [1 + {-1^5] S,,j ■ + aS,., • 5,+2,,} 

(17) 



2 J ' 



where Sij is the 3=1/2 spin operator at the ith site in 
the jth chain, J' is the interchain antiferromagnetic cou- 
pling, and denotes the summation over the pairs 
of nearest-neighbor chains. To put th roug h the inter- 
chain mean-field treatment [H IM IM IH Isl] . 
we introduce two kinds of mean fields induced by the in- 
terchain interaction and the external magnetic field as 
{Sf) = — (— l)*ms and (Sf) = m^. The former is the 
staggered magnetization and the latter is the uniform 
magnetization. The ID mean-field Hamiltonian thus ob- 
tained are described as 

Hmf - 2J^{[l + (-iy^]S,-5,;+i+aS,-S,+2} 

i 

- {gfiBH~hu)Y,S! ~hsY,S:, (18) 



where /in and hg are the effective internal fields given by 
hu = zJ'niu and hg = zJ'rUs-, respectively, with z being 
the number of the adjacent chains. The effective field is 
also defined as = gfisH — hu- We set J = gfj,B = 1. 



B. Staggered susceptibility 

On the basis of Hmf, we calculate the staggered mag- 
netizations per a site rus and the uniform magnetizations 
per a site by means of the infinite-system DMRG 
method. It is known that the original infinite-system 
DMRG method devised by White may suffer from a 
problem that the system in magnetic fields is trapped at 
metastable states during the renormalization process. To 
avoid this problem, we use the modified infinite-system 
DMRG algorithm by adopting the recursion relation to 
the wave function [43, . 

In given iJ^, Wg can be obtained as function of /ig, 
while in given z J', ttIu can be obtained as function of H^. 
The results are shown in Figs. 5 and 6. Note that for the 
pure ID system {,]' = 0), in a = 0.20 and 5 = 0.15 the 
longitudinal IC spin correlation becomes dominant (77^ < 
1 < 77^) around the half-magnetization of the saturation, 
while in a = 0.05 and 6 = 0.40 the transverse staggered 
spin correlation becomes dominant {rj^ < 1 < 77^) in 
< TO < 1/2. Solving the numerical results for iris self- 
consistently with hg = zJ'mg, we investigate whether the 
long-range staggered order exists or not. 
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FIG. 5: The staggered magnetization nis for several values of 
Hu as function of h^. Parameters used are (a) a = 0.20 and 
6 = 0.15, and (b) a = 0.05 and 5 = 0.40. 

We first discuss the results for a = 0.20 and S = 0.15 
shown in Fig. 5(a). At = 0, where the system has 



7 



an excitation gap, the magnetization curve of rris has a 
finite derivative at the origin, implying that the staggered 
susceptibility 



(19) 



takes a finite value. At = 4.60, niu shows the satura- 
tion as shown in Fig. 6 and the same behavior of Xs as 
in ifu = is observed. In 1.30 < < 3.99, the system 
becomes gapless. At = 1.30 and 3.99, Xs shows a di- 
vergence, indicating that the staggered long-range order 
is stabilized. In 1.70 < < 3.80, on the other hand, 
Xs takes finite values. The results implies that the stag- 
gered long-range order does not emerge up to a certain 
critical value zJ'^. For = 3.00, the critical value can 
be evaluated as zJ^ ^ 8.95 x 10~^, which is the largest 
in 1.70 < Hu < 3.80. The uniform magnetization rriu 
for a = 0.20 and 5 = 0.15 in Fig. 6 has been calculated 
using zJ' = 8.95 x 10~^. From the results for m^, we 
find that in 1.70 < < 3.80 the uniform magnetization 
varies in 0.10 < ruu < 0.32. 

In a = 0.05 and S = 0.40, the behavior of Xs is quite 
different as compared with that shown in Fig. 5(a). As 
shown in Fig. 5(b), Xs shows a divergence in 2.60 < 
Hu < 3.99, where the uniform magnetization varies in 
< rriu < 1/2. The results indicate that the staggered 
long-range order is always stabilized in < rriu < 1/2. 



method based on the Lanczos algorithm, we calculate 
the static structure factor. The results are shown in Fig. 
7(a). 



ct=0.20 
5=0.15 



• m=l/3 
^ m=l/4- 
■ m=l/6 





FIG. 6: The uniform magnetization niu for {a, 5) = 
(0.20, 0.15), (0.05, 0.40), and (0.15,0.43) as function of Hu- 
The interchain couplings are set to be zJ' = 8.95x 10~^ , 2.50x 
10"^, and 6.85 x 10"^, respectively. 



C. Static structure factor 

We investigate characteristics of the ordered state in 
1.70 < i?u < 3.80 for a = 0.20 and 5 = 0.15. Since the 
staggered long-range order does not appear in zJ < zJ^ 
in this case, the staggered magnetization becomes zero 
and then hg = 0. Using the numerical diagonalization 



10 20 30 40 

FIG. 7: (a) Static structure factor ^^(fc) in a = 0.20 and 
5 = 0.15 for ruu = 1/3, 1/4, and 1/6. (b) N dependence of 
S^(2A:f). The solid lines are obtained by the least-squares 
method. The divergence powers are evaluated as 0.196 for 
rUu = 1/6, 0.237 for niu = 1/4, and 0.170 for rUu = 1/3. 

From the magnetization curve for shown in Fig. 6, 
we estimate that = 1/3, 1/4, and 1/6 correspond to 
Fu - 3.88, 3.44, and 2.50, respectively. For these H^, the 
staggered long-range order docs not emerge in zJ < zJ^ 
as shown in Fig. 5(a). The size dependence oi S^{k) has 
an algebraic singularity at A: = 2kF with increasing the 
system size A/' as shown in Fig. 7(b). By the least-squares 
method, the power is evaluated as 0.196 for = 1/6, 
0.237 for mu = 1/4, and 0.170 for = 1/3. From 
the results, we conclude that the system has a tendency 
towards the formation of an IC long-range order in 1.70 < 
ifu < 3.80 for zJ < zJ^. 

We calculate rris, mu, and S^{k) in other parameter 
sets of a and S systematically using the phase diagram 
shown in Fig. 2. We find that in the darkly shaded area of 
Fig. 2 the system has a tendency towards the formation 
of an IC long-range order with the period l/rUu around 
TOu = 1/4 in zJ < zJ^. 
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IV. DISCUSSION 



FsPNN 



is adequate, Hp^ will be slightly larger than 3.77T. Of 
course, the parameters a and 6 adequate for l/Ti have 
to well reproduce the experimental results for thermody- 
namic quantities |49l l5l| . 



Recently, 1/Ti for F5PNN, which is considered to be 
a ID S* = 1/2 bond-alternating spin system was 
measured in magnetic fields 0]. In Hc^ = 2.5T < H < 
= 6.5T, 1/Ti exhibited power-law behavior. Fur- 
thermore, at H — 5.2T, l/Ti becomes almost indepen- 
dent of temperature. It seems difficult to explain the 
results for l/Ti at H = 5.2T on the basis of the S* = 1/2 
bond-alternating spin-chain model. We thus take ac- 
count of a small NNN interaction in addition to the bond 
alternation ISOl. 




FIG. 8: The divergence exponent 7 for F5PNN (q = 0.15 
and 5 = 0.43) as function of H, where l/Ti ~ T"^. Using 
the experimental data, Hp^ and Hp^ are evaluated as Hp^ ^ 
3.77T and Hp^ ~ 5.79T. The line is to guide the eyes. 

On the basis of the Hamiltonian we investigate 
the divergence exponent 7 of 1/Ti for F5PNN numeri- 
cally. Since the bond-alternation parameter of F5PNN 
was evaluated as S = 0.43 we slightly increase a in 
a fixed 6 — 0.43 and investigate whether the behavior 
r/^ < 1 < 77^ appears in < H < H^^- As shown 
in the dotted line in Fig. 2, such behavior emerges at 
a ^ 0.15. The field dependence of 7 is shown in Fig. 8. 
In iJci = 2.36 < H < = 4.00 the system is gapless, 
and at Hp-^ — 2.88 and Hp^ = 3.71, 1/Ti becomes inde- 
pendent of temperature. Using the experimental data of 
and , we evaluate the strength of Hp^ and Hp^ ; 
Hp, = 2.5T-K(2.88-2.36)/(4.00-2.36) x (6.5T-2.5T) - 
3.77T and Hp^ = 2.5T + (3.71 - 2.36)/(4.00 - 2.36) x 
(6.5T- 2.5T) - 5.79T. Therefore, 5.2T in F5PNN prob- 
ably corresponds to Hp^ . In order to improve quantita- 
tive agreement, a slightly smaller a than 0.15 may be 
adequate within S = 0.43. We have shown that a small 
NNN interaction may be essential in F5PNN. To con- 
firm this scenario, the magnetic field corresponding to 
Hp-, have to be observed, where 1/Ti becomes also inde- 
pendent of temperature. If a slightly smaller a than 0.15 
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FIG. 9: (a) The staggered magnetization rris for several values 
of -ffu as function of /is. Parameters are a = 0.15 and & — 0.43. 
(b) Static structure factor S^{k) for mu — 1/4. Inset: L 
dependence of 5'^(2fcp). From the solid lines obtained by the 
least-squares method, the divergence power is evaluated as 
0.346. 

For F5PNN, the transition into the 3D ordered state 
in magnetic fields was observed in further low tempera- 
tures by the specific heat measurements jHJ • We discuss 
the characteristics of the 3D ordered state in F5PNN on 
the basis of the Hamiltonian p8|l. The staggered suscep- 
tibility and the static structure factor are calculated for 
a — 0.15 and S ~ 0.43 by the same method as used in Sec. 
III. The results are shown in Fig. 9. In H^ = 3.30 and 
3.60, Xs takes finite values, implying that the staggered 
long-range order does not emerge up to a certain critical 
value zJ^. As shown in Fig. 6, the uniform magnetization 
takes the values niu = 0.22 for H^ — 3.30 and mu = 0.27 
for Hu — 3.60. Note that mu is calculated using the 
critical value 2 — 6.85 x 10~^ at H^ — 3.60, which is 
larger than that at H^ = 3.30. In the inset of Fig. 9(a), 
the size dependence of iS'^(2A:f) for mu — l/4(^ 0.22 and 
0.27) is shown. An algebraic singularity of S'^(2A:f) as 
L — > cx) is seen with the power 0.346. Therefore, the 
system has a tendency to form an IC long-range order 
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around rriu = 1/4 in zJ < zJ^. A detailed experimental 
study for the ordered state of F5PNN under magnetic 
fields is desirable. 

B. Dominant longitudinal IC spin correlations in 
other systems 

We have argued that a dominant longitudinal IC spin 
correlation is attributed to the formation of the half- 
magnetization plateau, which can be regarded as the 
CDW insulating state. Therefore, similar behavior is 
expected in other ID spin-gapped systems in magnetic 
fields, e.g., the 5 = 1/2 two-leg spin-ladder system with a 
cyclic four-spin interaction. In this system, it was shown 
that the CDW-like half-magnetization plateau appears 

For a S = 1/2 two- leg spin-ladder material SrCu2 03, 
two models were proposed to reproduce the tempera- 
ture dependence of the susceptibility. The one model de- 
scribes the pure two-leg ladder the other model 
includes the effects of a cyclic four-spin interaction [s^ . 
To distinguish characteristics between the two models, 
some experimental methods for the observation of dy- 
namical properties were proposed [s^ |5^ |53|- If the 
temperature-independent 1/Ti at certain two fields in 
< H < will be observed, such findings are also 
an evidence for the effects of a cyclic four-spin interaction 
in SrCu2 03. 

V. SUMMARY 

We have first investigated the critical properties of the 
5=1/2 bond-alternating spin chain with a NNN interac- 
tion in magnetic fields. From the results obtained by the 
numerical calculation and those obtained based on the 
effective Hamiltonian, we have concluded that there is a 
parameter region where the longitudinal IC spin corre- 



lation becomes dominant around the half-magnetization 
plateau. The results are regarded as a manifestation of 
the nature of the TL liquid close to the CDW transition 
point. Experimentally, such behavior can be observe by 
the field dependence of the divergence exponent of 1/Ti 
with decreasing temperature. 

When temperature is further decreased, the interchain 
interaction becomes relevant. We have calculated the 
staggered susceptibility and the uniform magnetization, 
combining the DMRG method with the interchain mean- 
field theory. In the parameter region where the dominant 
longitudinal IC spin correlation appears for J' — 0, the 
staggered long-range order does not emerge up to a cer- 
tain critical value zJ^ around ~ 1/4, while in the 
other parameter region, the staggered long-range order 
is stabilized in < TOu < 1/2. To investigate the char- 
acteristics of the long-range order in the former parame- 
ter region, we have calculated the static structure factor. 
From the size dependence of S^(2kp), we have shown that 
the system has a tendency to form an IC long-range or- 
der around ruu = 1/4 in 2; J < zJ^. Using the results, we 
have discussed the recent experimental results for 1/Ti 
in magnetic fields performed for F5PNN. 
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